Self-Duality in Superconductor- Insulator Quantum Phase Transitions 
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It is argued that close to a Coulomb interacting quantum critical point, the interaction between two vortices 
in a disordered superconducting thin film separated by a distance r changes from logarithmic in the mean-field 
region to 1/r in the region dominated by quantum critical fluctuations. This gives support to the charge- vortex 
duality picture of the observed reflection symmetry in the current-voltage characteristics on both sides of the 
transition. 
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One of the most intriguing results found in experiments 
on quantum phase transitions in superconducting films, 2- 
dimensional Josephson-junction aiTays [|l]], quantum Hall sys- 
tems 1^, and 2-dimensional electron systems ||^] is the strik- 
ing similarity in the current-voltage {I-V) characteristics on 
both sides of the transition. By interchanging the / and V 
axes in one phase, an I-V characteristic of that phase at a 
given value of the applied magnetic field (in superconduct- 
ing films, 2-dimensional Josephson-junction arrays, and quan- 
tum Hall systems) or charge-caiTier density (in 2-dimensional 
electron systems) can be mapped onto an I-V characteristic 
of the other phase at a different value of the magnetic field 
or charge-carrier density. This reflection symmetry hints at 
a deep connection between the conduction mechanisms in the 
two phases that can be understood by invoking a duality trans- 
formation [0^]. Whereas the conducting phase is most suc- 
cinctly described in terms of charge carriers of the system, the 
insulating phase is best formulated in terms of vortices, which 
behave as quantum point particles in these systems. The du- 
ality transformation links the two surprisingly similar looking 
descriptions. 

There appears to be, however, one disturbing difference. 
Whereas charges interact via the usual 3 -dimensional 1/r 
Coulomb potential, vortices are believed to interact via a log- 
arithmic potential-at least for distances smaller than the trans- 
verse magnetic penetration depth A^, which is typically larger 
than the sample size [|[]. This is disturbing because the differ- 
ence should spoil the experimentally observed reflection sym- 
metry. 

It is this fundamental problem we wish to address in this 
Letter. It will be shown that close to a Coulomb interacting 
quantum critical point (CQCP), the interaction between vor- 
tices in disordered superconducting films changes from log- 
arithmic in the mean-field region to 1 /r in the region domi- 
nated by quantum critical fluctuations. This conclusion is an 
exact result, depending only on the presence of a CQCP. 

A common characteristic of the systems mentioned is, apart 
from impurities, the presence of charge carriers confined to 
move in a 2-dimensional plane. As the 1/r Coulomb repul- 
sion between charges is genuine 3-dimensional, we assume 
this interaction not to be affected by what happens in the film, 
which constitutes a mere slice of 3-dimensional space. In con- 



trast to this, the interaction between vortices is susceptible to 
the presence of a CQCP. This is because the vortex interac- 
tion is a result of currents around the vortex cores which are 
confined to the plane. 

As starting point, we take the observation (for a review, 
see Ref. [^) that close to a CQCP, the electric field E scales 
with the correlation length ^ as £' Cr^C~^ ~ £,~^^^^\ 
Here, denotes the correlation time, indicating the time pe- 
riod over which the system fluctuates coherently, and z is the 
dynamic exponent. Thus conductivity measurements [||,^ 
close to a CQCP collapse onto a single curve when plot- 
ted as function of the dimensionless combination (5''(^+i) /E, 
where 5 — {K — K^jK^ measures the distance from the 
critical point K^, and v is the correlation length exponent, 
^ ^ . (For a field-controlled transition, K stands for the 
applied magnetic field, while for a density-controlled transi- 
tion it stands for the charge-carrier density.) The scaling of 
the electric field with the correlation length expresses the more 
fundamental result that the anomalous scaling dimension cJa 
of the magnetic vector potential A is unity, (Za = 1- 

In addition, because the magnetic vector potential always 
appears in the gauge-invariant combination V — gA, the 
anomalous scaling dimension of the electric charge q of the 
charge carriers times the vector potential is unity too, d^A = 
1. Writing the anomalous scaling dimension of the vector po- 
tential as a sum ^a = c^a + ^^?a of its canonical scaling 
dimension + z — 2), obtained by simple power 

counting, and (half) the critical exponent 77A, describing how 
the correlation function decays at the critical point, we con- 
clude that dq = d*^ — ^TjA- Here, ~ 1 ~ '^a stands for 
the canonical scaling dimension of the electric charge. Now, 
for a 1/r Coulomb potential, the charge scales as ~ 
independent of the number d of space dimensions [^] . Com- 
bined with the previous result, this fixes the value of the criti- 
cal point decay exponent r/A in terms of the number of space 
dimensions and the dynamic exponent: 



77A = 5 — d — 2z. 



(1) 



In Ref. it was further argued that in the presence of dis- 
order, the electric charge is finite at a CQCP, so that z = 1. 
This prediction was first confirmed for disordered supercon- 
ducting films []lO[], and subsequently also for 2-dimensional 
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Josephson-junction arrays [|l]], quantum Hall systems [[IJ 
and 2-dimensional electron systems [^. With z = 1, the 
value of the critical-point decay exponent becomes tja = 1 
in d = 2. As we will now demonstrate, this leads to a quali- 
tative change in the interaction potential between two vortices 
from logarithmic in the mean-field region, where t^a = 0, to 
l/r in the vicinity of the CQCP, where t^a = 1- 

To set the stage, let us first consider a bulk superconductor 
with two static vortices directed along the x^-axis and sepa- 
rated a distance r. For our purposes, the effective phase-only 
I p^ ] Hamiltonian HeS = (Ps/2m^)(V(^ - gA)^ in terms 
of the phase (p of the superconducting order parameter — the 
so-called Anderson-Bogoliubov mode — suffices (for reviews, 
see Ref. [^). Here, ps is the superconducting mass density, 
which scales as ps ^ ^^-(d+z) [[j^^ ^jj^j jg jjjg mass of 
the charge carriers. The interaction potential can be extracted 
from the magnetic part of the effective action S'mag ■ Written 
as a functional integral over the the magnetic vector potential, 
it is given in the Coulomb gauge V • A = by 



(2a) 



with A the magnetic penetration depth, which is related to ps 
via A^^ = q^ps/m'^. The mass term is generated through 
the Anderson-Higgs mechanism by integrating out the phase 
mode (fi. The so-called plastic field [Q 
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(2b) 



with $0 ~ 27r/g the magnetic flux quantum in units where the 
speed of light and Planck's constant h is set to unity, describes 
the two vortices located along the lines Ca (a — 1, 2). 

Note that since the anomalous scaling dimension of the 
magnetic vector potential is unity, the dimension of the 
Maxwell term is 4, implying that in d = 2 it is an irrele- 
vant operator in the renormalization-group sense. This term 
is, however, important when considering the interaction be- 
tween vortices. 

To facilitate the calculation in the case of a superconducting 
film below, we linearize the first term in Eq. ( Pa| ) by introduc- 
ing an auxiliary field h via a Hubbard-Stratonovich transfor- 
mation to obtain the combination i(V x A — B^) • h — |h^. 
After integrating out the magnetic vector potential, we arrive 
at a form appropriate for a dual description in terms of mag- 
netic vortices rather than electric charges [ [l6| ] 



^(Vxh)2-ih"-ihBP 



(3) 



Physically, h represents (i times) the fluctuating local induc- 
tion; it satisfies the condition V • h = 0. The vortices couple 
with a coupling constant g = $0/ A independent of the electric 
charge to h. Observe the close similarity between the origi- 
nal ( pa| ) and the dual form (^j). This becomes even more so 
when an external electric current is coupled to the A field 
by including a term —A ■ in Eq. (^, and B^ describing 
the vortices is set to zero there. 



Integrating out the local induction, one obtains the well- 
known Biot-Savart law for the interaction potential 5mag — 
dtV between two static vortices in a bulk superconductor 



Vir) = ^J d3xd3yBf(x)G(x-y)i?f(y) 
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-|-i[ln(r/2A)+7]+0(r/A)2 



(4) 



where we ignored the self-interaction. In Eq. (^, G'(x) is the 
correlation function whose Fourier transform reads G(k) = 
l/(k^ + A^^), R denotes the distance between the differen- 
tial lengths dl^ and dP, L is the length of each of the two 
vortices, and 7 is Euler's constant. For distances smaller than 
the magnetic penetration depth, which is the length scale for 
variations in the current and the magnetic field, the interaction 
is logarithmic as in a superfluid. If the system size is smaller 
than A, it will replace A as infra-red cutoff in the logarithm, 
and there will be no reference to the electric charge anymore. 

To describe magnetic vortices in a film of thickness w JT^], 
the bulk result (^ has to be adjusted in two ways to account 
for the fact that both the vortices and the screening currents, 
which produce the second term in (^, are confined to the 
plane. This is achieved by including a Dirac delta function 
'w5{x'i) in the second and third term. Instead of Eq. (^, one 
then arrives at the interaction potential [ |l8| , p^ 



27r 



[ln(^/4Ax)+7]+0('VAi)', 



(5a) 



where B^ = — $0 '^(^^ ^ ^1) describes the vortices in 
the film with coordinates x^, A^ = X^/w is the transverse 
magnetic penetration depth, g'j^ = $q / A^ the coupling con- 
stant squared, and 



G±{x±) = J dx3G±{x±,x3) 



d^fcx 
(2^ 



-G±(ki,0), 



(5b) 



with G_L(kx,0) = 2/fcx(2fc^ + AJ^). For smaU distances, 
the interaction is seen to be identical to that in a bulk super- 
conductor [p^, and also to that in a superfluid film. As in the 
bulk, the vortex coupling constant g± in the film is indepen- 
dent of the electric charge, g'j^ = <I>g/A^ = {2^)^ psW / m? , 
with Ps the bulk superconducting mass density. 

The above results are valid in the mean-field region, where 
77A = 0. In the critical region governed by a CQCP, the value 
of this exponent is unity, and the correlation function becomes 



Gi(ki,0) = 



(6) 
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with ^ the field renormalization factor. Because the 
magnetic vector potential and the local induction renormalize 
in the same way, their renormalization factor is identical. Due 
to this extra factor, the interaction between two vortices in the 
film takes the form of a 1/r Coulomb potential 



gia 

27r r ' 



(7) 



where a is some microscopic length scale which accompanies 
the renormalization factor Za for dimensional reasons 1 19 1. 



Since the electric charge is finite at the CQCP, the pene- 
tration depth X± oc 1/ps scales with the correlation length 
as ^ In the correlation function (^) we thus have the 
combination l/{2k± + ^^^) which should be compared with 
l/(k^ + ^~^) for a bulk superconductor 

The absence of any reference to the electric charge in the 
renormalized and bare interaction (at least for small enough 
systems) implies that the same results should be derivable 
from our starting Hamiltonian restricted to two dimensions 
and with g set to zero: H± = {psw/2m'^){V ±ip — (f^)'^. The 
plastic field ip^, with x (p^ = — 27r^^(5(x^ — x") 
describes vortices in a superfluid [|l5]]. It is obtained from 
the description involving the plastic field by a canonical 
transformation of the vector potential. By directly integrat- 
ing out the Anderson-Bogoliubov mode, and ignoring the k± 
dependence of ps, which is valid outside of the critical re- 
gion, one easily reproduces the bare interaction potential (^). 
The renormalized interaction (Q) is obtained by realizing that 
the anomalous scaling dimension of the superconducting mass 
density is = {d + z) — 2 [p4|], so that in our case p^ ^ kj_. 
In other words, the extra factor of k± that came in via the 
renormalization factor Za in our first calculation to produce 
the 1 /r potential, comes in via ps here [^0|]. 

A similar change in the r-dependence of the interaction be- 
tween two vortices upon entering a critical region has been 
observed numerically in the 3-dimensional Ginzburg-Landau 
model [gip. Near the charged fixed point of that theory, 
?7a = 1 [p2[], as in our case. 

This is a very pleasing coincidence as the (2+1)- 
dimensional Ginzburg-Landau model constitutes the dual for- 
mulation of the system. To appreciate the basic elements of 
the dual theory, note that the dynamics of the charged degrees 
of freedom is described by the effective Lagrangian 
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2to2 



(8) 



with c the speed of sound. In accord with the above findings, 
we have ignored the coupling to the magnetic vector potential, 
so that Eq. ^ essentially describes a superfluid. Although the 
complete effective theory is Galilei invariant | p3| , p4[ ], the lin- 
earized form is invariant under Lorentz transformations, 
with c replacing the speed of light. 

In the dual formulation, where the roles of charges and vor- 
tices are interchanged, the Anderson-Bogoliubov mode me- 
diating the interaction between two vortices is represented as 
a photon associated with a fictitious gauge field a^, i.e., (in 
relativistic notation) d^i^ ^ €^^\d'^a^. In 2h-1 dimensions. 



a photon has only one transverse direction and thus only one 
degree of freedom — as has the Anderson-Bogoliubov mode. 
The elementary excitations of the dual theory are the vortices, 
described by a complex scalar field ip. Specifically, the (well- 
known) dual theory of Eq. (ph is the Ginzburg-Landau model 



WW 

(9) 



with /^^ — d^a^ — d^a^, a mass parameter, and u the 
strength of the self-coupling. Both the gauge part as well as 
the matter part of the dual theory are of a relativistic form. 
The gauge part is because the effective theory ^ is Lorentz 
invariant, while the matter part is because vortices of positive 
and negative circulation can annihilate, and can also be cre- 
ated. In this sense they behave as relativistic particles. As was 
pointed out in Ref. [@], the speed of "light" in the gauge and 
matter part in general differ. 

The interaction potential (|^ between two external vortices 
is now being interpreted as the 2-dimensional Coulomb po- 
tential between charges. The observation concerning the crit- 
ical behavior of the Ginzburg-Landau model implies that the 
qualitative change in V{r) upon entering the critical region is 
properly represented in the dual formulation. 

Whereas in the conducting phase, the charges are con- 
densed, in the insulating phase, the vortices are condensed 
[Q]. In the dual theory, the vortex condensate is represented 
by a nonzero expectation value of the ?/; field, which in turn 
leads via the Anderson-Higgs mechanism to a mass term for 
the gauge field a^. Because {e^i^xd'^ Y ^ (9^(^)2, the mass 
term with two derivatives less implies that the Anderson- 
Bogoliubov mode has acquired an energy gap. That is to say, 
the phase where the vortices are condensed is indeed an in- 
sulator. Since electric charges are seen by the dual theory as 
flux quanta, they are expelled from the system as long as the 
dual theory is in the Meissner state. Above the critical field 
h = X a = hc^ they start penetrating the system and 
form an Abrikosov lattice. In the original formulation, this 
corresponds to a Wigner crystal of the charges. Finally, when 
more charges are added and the dual field reaches the critical 
value /ic2, the lattice melts and the charges condense in the 
superfluid phase described by the effective theory dq). 
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